A Shallow Water Analogue of the Standing Accretion Shock Instability: 
Experimental Demonstration and Two-Dimensional Model 
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Despite the sphericity of the collapsing stellar core, the birth conditions of neutron stars can be 
highly non spherical due to a hydrodynamical instability of the shocked accretion flow. Here we 
report the first laboratory experiment of a shallow water analogue, based on the physics of hydraulic 
jumps. Both the experiment and its shallow water modeling demonstrate a robust linear instability 
and nonlinear properties of symmetry breaking, in a system which is one million times smaller and 
about hundred times slower than its astrophysical analogue. 
PACS numbers: 97.60. Bw, 95.30.Lz, 47.20.-k, 01.50.Pa 



The physics of stellar core-collapse involves complex in- 
gredients such as nuclear physics, neutrino interactions, 
multidimensional hydrodynamics and general relativity. 
Current understanding is based on simplified formula- 
tions which make the problem numerically tractable at 
least, if not physically intuitive. A recent breakthrough 
occurred on the hydrodynamical side with the recogni- 
tion of a new instability mechanism [H-0] taking place 
in the inner 200km of the collapsing core, during a few 
hundred milliseconds before the explosion, while a spher- 
ical accretion shock stalls above the surface of the proto- 
neutron star. This shock was found to be unstable with 
respect to global oscillations, with a period of a few tens 
of milliseconds. An explosion is launched when sufficient 
neutrino energy diffusing out of the dense central region 
has been intercepted by the post-shock accreting mat- 
ter When the explosion proceeds, the large scale 
asymmetry resulting from this standing accretion shock 
instability (SASI) imposes on the neutron star a velocity 
kick mediated by gravitational forces [Ml , and can also 
affect its spin in a direction surprisingly opposite to the 
global rotation of the shock wave [IH H] • The pulsar spin 
period that could be accounted for by the spiral mode of 
SASI is still debated |12Hl5l |. These effects are likely to 
depend on the progenitor mass and rotation rate, the 
physical ingredients of the numerical model, and could 




also be sensitive to numerical limitations such as numer- 
ical resolution. We shed new light on the physical nature 
of these processes by adopting an experimental approach 
for the first time. The outcome is strikingly similar to the 
astrophysical results obtained through numerical simula- 
tions in an adiabatic approximation [l3j |. The simplicity 
of the experiment in the familiar environment of shallow 
water physics should allow physicists to build up physical 
intuition about this instability. 

A Shallow Water Analogue to a Shock Instability. The 
experimental set up named SWASI is designed to produce 
a Shallow Water Analogue of a Shock Instability. It is 
based on the analogy between acoustic waves in a gas and 



FIG. 1. Like the classical hydraulic jump in a kitchen sink 
(upper left), the SWASI experiment (upper right and bot- 
tom) involves a hydraulic jump associated to the deceleration 
of a radial flow of water. Water is injected inward from an 
annular injection reservoir (R) along a hyperbolic potential 
well, and evacuated through a vertical cylinder (C), whose 
walls mimic the surface of the neutron star. A pump (P) 
distributes collected water. The lower picture illustrates the 
stationary and axisymmetric character of water injection in a 
stable configuration. 



surface gravity waves in shallow water [16| . Their veloc- 
ity c when the depth H is smaller than their wavelength 
is defined by c 2 = gH, where g is the gravitational accel- 
eration in the laboratory. The hyperbolic shape H greLV (r) 
of the surface on which the shallow fluid flows is a for- 
mal analogue to the gravitational potential $ = g-ff grav , 
dominated by the mass of the central proto- neutron star. 
The system of equations defining this flow in the inviscid 
2D approximation can be directly compared to the equa- 
tions describing an isentropic gas with an adiabatic index 
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7 = 2 [13], where the compressible gas density plays the 
same role as the depth of the fluid layer in the experi- 
ment: 
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The Froude number Fr = \v\/c comparing the fluid ve- 
locity v to the wave velocity plays the same role as the 
Mach number in a g as. Shock waves are analogue to 
hydraulic jumps 16j. Both can be idealized as sharp 
discontinuities where the mass flux and momentum flux 
are conserved. In contrast with the kitchen sink experi- 
ment |18| (Fig. [TJ upper left), the fluid is injected inward 



in the SWASI experiment (Fig. [TJ bottom). The inner 
boundary is made of a vertical pipe such that water ac- 
cumulating in the gravitational potential is continuously 
extracted by spilling over its upper edge (Fig. [TJ upper 
right). The height of the upper edge of the inner cylinder 
defines a pressure threshold over which water is efficiently 
evacuated. It is adjusted vertically in order to choose the 
radius of the hydraulic jump (rj p ~ 20cm). Conveniently 
simple, this inner boundary could be interpreted as an 
analogue of a cooling process such as neutrino emission 
close to the neutrinosphere, which decreases the energy 
of the gas when the density and temperature are high 
enough. The experiment size was chosen large enough 
in order to minimize the effects of viscosity: the viscous 
drag is negligible except in the thin layer ahead of the 
jump. The injection radius is 32cm, and the radius of 
the inner cylinder is = 4cm. The shape of the gravita- 
tional potential is described by di? grav /dr = (5.6cm/r) 2 . 
Accreted water is pumped back into the annular reservoir 
through 16 pipes. A layer of sand dampens any inhomo- 
geneity before injection. 

Despite its simple set up, this experiment is expected 
to capture some hydrodynamical properties of the accret- 
ing gas in the equatorial plane of the stellar core 3, 3] , 
particularly those observed in the adiabatic approxima- 
tion [13J. The shallow water model can be scaled using 
the jump radius rj p as distance unit, the free fall veloc- 
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ity Vfi = (2gH 3 ^ av ) 1/2 as velocity unit, and the free fall 
timescale 



fjp/uff as time unit. The dimensionless 
solution depends on three parameters only [17] : the rel- 
ative size of the inner boundary r^/r^-p, the pre-shock 
velocity vi/vg and the Froude number Fri ahead of the 
shock. The experimental results can be scaled to astro- 
physical proportions by using the ratio r s h/rj p ~ 10 6 for 
distances, and the ratio tg/tg for timescales: 
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where Mns ~ 1.2M so i is the mass of the proto-neutron 
star and G is the gravitational constant. 



Surface gravity waves and advected vorticity pertur- 
bations present in shallow water are directly compara- 
ble to acoustic and vorticity waves in a compressible gas. 
These are a possible source of an unstable cycle similar to 
the advective-acoustic cycle in SASI 0-0| and also com- 
parable to the vortical-acoustic cycle seen in a shocked 
isothermal gas [20j |. 

Symmetry breaking. Indeed, as the water flux is in- 
creased, a large scale instability sets in through growing 
oscillations of the hydraulic jump (Fig. [21 left). This in- 
stability is also observed in numerical simulations 21] 
of the experimental setup in the 2D shallow water ap- 
proximation (Fig. [21 right). The visual resemblance with 
astrophysical simulations is supported theoretically by 
the formal similarity between the set of shallow water 
equations and the set of adiabatic gas equations used in 
13]. A perturbative analysis of the shallow water equa- 
tions [17[ reveals a SASI-like instability dominated by the 
global mode m=l in most of the parameter space (Fig. [4]). 
The oscillation period measured in the experiment is in 
excellent agreement with both the perturbative analysis 
and the numerical simulations (Fig. [3]). Using the scal- 
ing factor in Eq. ([3]), an oscillation frequency of 3s and 
a growth rate of 0.2s _1 in the experiment for v\/vs ~ 1, 
r jp/ r * = ^ an d F r i ~ 4 would correspond to an oscilla- 
tion period of 42ms and a growth time of 70ms on an 
astrophysical scale, which are comparable to the values 
measured in core-collapse models (e.g. Fig. 17 in Q). 

Like 3D astrophysical simulations [l3, 14], the exper- 
iment displays both the sloshing and the spiral modes 
(movies Ml, M2, M3 [2l|). These are easily triggered de- 
pending on the random initial perturbations, as expected 
since they have the same linear growth rates. However 
they interact during the nonlinear phase of their evolu- 
tion, ultimately favoring a single right or left spiral mode 
even though the injected flow contains no angular mo- 
mentum. This nonlinear behavior, first noted in astro- 
physical numerical simulations [13], is clearly observed 
in the experiment if the instability is vigorous (movies 
M2, M3 and M5 2l|). The shape and dynamical evo- 
lution of the hydraulic jump in the nonlinear regime is 
remarkably akin to the astrophysical numerical simula- 
tions HI, 3, 0, EI , with the formation of a triple point 
when the spiral mode reaches nonlinear amplitudes (Fig. 
[21 movie M6 [2l|). Another striking similarity is the fact 
that the angular momentum of the accreted fluid is op- 
posite to the direction of rotation of the hydraulic jump. 
This is made visible in the experiment by using a light 
horizontal bar freely rotating on a vertical axis, with both 
ends bent vertically and immersed in the flowing water. 
The hydraulic jump and the horizontal bar systematically 
rotate in opposite directions (Fig. [2] and movies M4, M5 

[HI). 

Effect of the bottom drag for a viscous fluid. The 
Reynolds number in the experiment can be defined using 
the viscosity v of water, the fluid velocity v and depth 
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FIG. 2. The hydraulic jump governed by a spiral mode dis- 
plays a rotating triple point both in the experiment (left) 
and in the numerical simulation of 2D shallow water equa- 
tions (right). The altitude of the free surface (upper right) 
and the vorticity (lower right) are shown. The same shape 
and dynamics are observed in astrophysical simulations 
The angular momentum in the accreted flow is visualized by 
a horizontal bar in the experiment (lower left). It spins in 
the direction opposite to the hydraulic jump (movie M4 [2l]|). 
The vorticity trails shown in the numerical simulation (lower 
right) of the experiment illustrate these counterrotating mo- 
tions (movie M6 
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FIG. 3. The oscillation period measured in the experiment, 
shown here as a function of the jump radius, is well cap- 
tured by both the perturbative analysis and the numerical 
simulation. Error bars relative to the measurement of the 
jump radius (±5mm) and the oscillation period (±3.5%) are 
indicated. The measured flow rate is lL/s (±5%) and the 
injection slit size is 0.74mm (±0.02mm). 
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FIG. 4. The stability of the flow depends on the strength of 
the hydraulic jump measured by the pre-jump Froude num- 
ber Fri , and the pre-jump velocity v\ measured in units of the 
free-fall velocity vg along the potential surface. The growth 
rate of the dominant m = 1 mode, deduced from the pertur- 
bative analysis of the 2D shallow water model, is displayed 
by contour lines for 0, 0.035, 0.07, 0.105, 0.14, 0.175 (f^) -1 , 
for a jump radius r- JP /r* — 4.9 (tg = 0.35s). The symbols '+' 
indicate the set of parameters used in experiments showing 
a robust m = 1 instability, the filled square corresponds to 
the parameters of Fig. [3] Dashed lines of constant flow rate 
(0.5, 1.5 and 2.5 L/s) from bottom to top correspond to a 
Reynolds number 408, 1224 and 2040 respectively, measured 
at the jump radius. 



H. It is smallest in the outer part of the experiment 
and inversely scales with radius: IZe = Hv/v — Q/2ni>r 
where Q is the injected flow rate. According to Fig. [H the 
Reynolds number at the shock radius is in the range 500- 
2000 within the parameters of the experiment, which sup- 
ports the laminar approximation for the viscous drag [23j , 
proportional to v/H 2 . This viscous drag is non negligi- 
ble only in the most shallow layer ahead of the hydraulic 
jump. The strength of the viscous drag has been esti- 
mated experimentally in a stationary supercritical flow 
(Fr> 1) by measuring the radial profile of the potential 
surface and the free surface of the fluid between R = 6cm 
and R = 31cm for an injection slit size of 1.3mm and a 
flow rate of 1.15L/s. The corresponding velocity pro- 
file was best fitted with an effective viscosity coefficient 
v = 0.03 cm 2 /s, such that the viscous drag is written as 
an additional term vv /H 2 in Eq. ©. The perturbative 
analysis [17| of the 2D shallow water model assumes that 
this parameter is a constant. 

In order to stress the analogy between the dynamics of 
the post-jump region in the experiment and the dynamics 
of the post-shock region in the astrophysical simulations, 
the flow parameters were chosen according to their values 
at the shock (velocity and Froude number) rather than 
at injection. Measuring the pre-jump velocity in units of 
the inviscid free fall velocity in Fig. [U the contour lines 
showing the growth rate of the instability for v — 0.03 
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cm 2 /s are well approximated by those obtained for an 
inviscid flow. 

Modeling approximations. The hydraulic jump has 
been modeled as a simple discontinuity where energy 
is dissipated through a viscous roller of negligible radial 
extension. In reality, depending on the strength of the 
jump, small scale waves can be generated which may in- 
terfere with the global instability for some parameters. 
Intermittent behavior of hydraulic jumps are well known 
in the range 2.5 < Fr < 4.5, while steadier jumps are 
expected in the range 4.5 < Fr < 9 24]. This latter is 
more typical of pre-shock Mach numbers relevant to core- 
collapse studies. The experiments showed possible inter- 
ferences between the intermittent character of the jump 
and the global instability for Fr < 3. Fig. 0] shows the 
set of most robust experiments performed for Fr > 3.3. 
The transition to turbulence, expected around Re= 2000, 
would require a different prescription for the viscous drag 
(proportional to v 2 / H). Viscosity in the horizontal plane 
and surface tension have been neglected. The vertical 
structure of the flow has been idealized in order to ob- 
tain a 2D description of the experiment. A more accurate 
modeling of this vertical structure could be obtained by 
introduc ing the Coriolis and Boussinesq correction coef- 
ficients 25( for the kinetic energy and the momentum 
equations, in the St Venant system and in the bound- 
ary conditions. Such quantitative improvement would 
be at the expense of the simplicity of the analogy with 
gas equations. 

A new tool for supernova physics. Drawing a paral- 
lel between stellar explosions and shallow water physics 
is unexpected, but is deeply rooted in the universality 
of the laws of fluid mechanics. The complementarity 
between the experimental and numerical approaches is 
beneficial because their inherent limitations are differ- 
ent. The experiment together with its idealized modeling 
demonstrated the robustness of the SASI instability in a 
converging flow, even in a fluid without entropy gradi- 
ents. This emphasizes implicitly the importance of vor- 
ticity perturbations, like in the vortical-acoustic cycle in 
an isothermal flow [20j. This new tool will be used to 
characterize nonlinear effects determining both the satu- 
ration amplitude and the interaction between right and 
left spiral modes. The saturation amplitude in the invis- 
cid approximation is expected to be governed by the par- 
asitic growth of the Kelvin- Helmholtz instability [26| , ap- 
parent in the numerical simulations (Fig. [2] lower right). 
Furthermore, a global rotation of the experiment will be 
introduced in order to mimic stellar spin. This is ex- 
pected to accelerate the growth of the spiral mode rotat- 
ing in the same direction as the progenitor [l^, 27 1 and 



temperature and pressure, and low cost. It is both a new 
research instrument for supernova science and an original 
tool for public outreach. 

The design and building of the experiment benefited 
from the help of P. Mulet, E. Gregoire, L. Rodriguez, F. 
Boissonnet and V. Foglizzo. This work is part of ANR 
funded projects Vortexplosion ANR-06-JCJC-0119 and 
SN2NS ANR-10-BLAN-0503. J.G. acknowledges support 
from STFC. 



enhance the influence of SASI on the angular momentum 
budget of the accretor. 

This shallow water analogue of a shock instability be- 
longs to the rare category of physical processes which are 
astrophysically relevant, accessible on Earth at normal 
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